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We analyze the consequences resulting from the asymmetric boundary conditions 
imposed by a non-uniform external magnetic field at the extremities of a planar 
Josephson tunnel junction and predict a number of testable signatures. When the 
junction length L is smaller than its Josephson penetration depth Aj, static ana- 
lytical calculations lead to a Fresnel-like magnetic diffraction pattern, rather than 
a Fraunhofer-like one typical of a uniform field. Numerical simulations allow to in- 
vestigate intermediate length (L ^ \j) and long [L > Xj) junctions. We consider 
both uniform and (5-shaped bias distributions. We also speculate on the possibility 
of exploiting the unique static properties of this system for basic experiments and 
devices. 

PACS numbers: 03.70.+k, 05.70.Fh, 03.65.Yz 

INTRODUCTION 

Both the static and dynamics properties of a Josephson tunnel junction (JTJ) are af- 
fected by the presence of an externally applied magnetic field. Since the discovery of the 
Josephson effect in 1962, the magnetic diffraction phenomena of the supercurrent and the 
occurrence of current singularities (Fiske steps) due to resonant cavity modes have been 
studied under the assumption of an homogeneous magnetic field [Ij. The aim of this letter 
is to overcome this textbook assumption and to investigate the consequences of enforcing 
asymmetric boundary conditions (ABC) on rectangular JTJs having length L along the X 
axis, width W along the Y axis and uniform critical current density Jc- We assume the 
junctions to be one-dimensional, i.e., W « L,Xj, with Xj — y^^^JT^Tr^e/Wc being the 
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Josephson penetration depthfS] (where $o is the magnetic flux quantum, de the effective 
junction magnetic penetration j3] and fio the vacuum permeability). It is well known that, 
under these circumstances, the gauge-invariant phase difference of the order parameters 
of the superconductors on each side of the tunnel barrier obeys the static or d.c. perturbed 
sine-Gordon equation[lJ: 

-0W = A|^ + ^, (1) 
in which the term Ib{X) is the distribution of the externally applied bias current 1^, choosing 

T /9 

the X-axis origin in the center of the junction, — J_^^^Xi){X)dX . In normalized units of 
X — X/Xj, Eq.([T]) becomes: 

02,^ = sin0(x) - 7(x). (2) 

With such normalization the junction length is / = L/Xj; further, the Josephson (zero- 
voltage) current ij through the barrier is obtained as the spatial average of sin0, ij — 

J /9 

(sin0) = (1/0 /_;^2 ^^^^(^)'^^' while the total bias current it — Ib/{JcWL) is given by the 
spatial average of 7, ib — (7). In this work we will consider two quite different symmetric 
bias current profiles: i) uniform bias j{x) — 7^ for which — and ii) (5-shaped bias 
7(3;) = ^sS{x) for which, being J"^ S{x)dx — 1, we have ib — ^s/l. (5-biased JTJs in a uniform 
magnetic field have been the subject of a recent theoretical and experimental investigation jH 
[5]. According to the magnetic Josephson equation[2j, the boundary conditions for Eq.([2]) 
depend on the external field values at the junction extremities: 

0,(-//2) = hL 0.(//2) = hR, (3) 

in which hL,R are the K-components Hlr of the applied magnetic field at the left and 
right junction ends, normalized to ^0/ (27iiiodeXj). Classically one considers the case of a 
uniform field applied in the junction plane perpendicular to the long junction dimension L, 
so that — Kl- We will refer to the above conditions as symmetric or classical boundary 
conditions to distinguish them from the asymmetric boundary conditions (ABC) achieved 
when the magnetic fields at the junction extremities have the same amplitudes, but opposite 
directions, that is, Kr — —hi = ha- We will analyze the properties of ABC later on; for 
now we only remark that they reminds of the boundary conditions for an in-line symmetric 
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j unction [0] with Kl^r being the self-fields produced by the external current which enters 
the junction at one extremity and leaves at the opposite one. Figs. 1(a) and (b) show two 
examples of how the ABC could be realized in practical devices. In both cases the asymmetry 
is obtained by the current Id flowing in a properly designed control line. The control line 
technique has been widely and successfully used to produce local magnetic fields in Josephson 
structures since 1994[7]. In Fig. 1(a) the control line, separated from the junction by an 
insulating layer, runs aside the long dimension of an overlap-type planar JTJ and fiips sides 
in the center of the junction. In the second, less easily achievable, case, sketched in Fig. 1(b), 
the control line is perpendicular to the junction plane and goes through a small hole drilled 
in the substrate in between the gap of a ring shaped JTJ. In the former case the bias current 
density is uniform, while in the latter case the current distribution can be expressed in terms 
of a (5-function[4j. It is quite evident that ABC are achieved in both electrode configurations. 




(a) (b) 

FIG. 1: Sketches of planar Josephson tunnel junctions in the asymmetric magnetic field Hd gen- 
erated by properly designed control lines, (a) (uniformly biased) linear and (b) ((5-biased) gapped 
annular Josephson tunnel junctions. The base electrode is in dark gray, the top electrode is in light 
gray and the junction area is white. 



ASYMMETRIC BOUNDARY CONDITIONS 

The ABC have several important implications. We begin by computing the Josephson 
current ij carried by the junction. From Eq.([2|) we easily get: 
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ij = y [0.(//2) - U-l/2)] + = ^'^^y^ + I,. (4) 

With the classical boundary conditions, the first term of the right side of Eq.Q vanishes, 
so that ij — ib- However, with ABC we get: 

= — + lb, (5) 

i.e., the zero- voltage current also depends on the amplitude of the asymmetric magnetic field. 
More precisely, the applied magnetic field acts as an extra field-dependent d.c. current source 
ih — '2ha/l in parallel to the external current source ib- In real units the normalized current 
ih corresponds to a supercurrent Ih — ihJc^L — '^HrW proportional to the junction width 
W (e.g., with Hji — lOA/m and W — 3/im, then Ih — GOfiA). Eq.Q suggests that a planar 
JTJ can be used as a magnetic first order gradiometer based on the readout of the current 
ih supplied by a stand-alone (unbiased) JTJ. However, it should be realized that, since ih is 
a zero- voltage current, the internal impedance of the current source is null, so that no power 
can be delivered to an external resistive load. 

From a mathematical point of view, the peculiarity of the ABC is that, provided that 
the current distribution is symmetric 7(— a:) = 7(a:), the solutions of Eq.([2]) have to be 
even functions, 0(x) = (j){—x), meaning that the phase difference between the two junction 
extremities is always null: 4>{l/2) — 0(— //2) = 0. As a consequence of the phase parity, we 
could limit our analysis to the spatial range [0,//2], with the condition in the origin to be 
found as follows. From Eq.(j2]), for any < < //2, we can write: 

/XO rxo 
sm(l){x)dx — / '){x)dx. (6) 
-XO J -XO 

Taking the limit Xq 0, the first integral vanishes, since 0(x) is a continuous function. If also 
the bias current density ^{x) is continuous, then the second integral vanishes too, enforcing 
0a;(O+) — 0x(O— ) = 0. Being the phase derivative an odd function, 0a;(O+) + 0a;(O— ) = 0, it 
vanishes in the origin: 



0x(O) = 0. (7) 

However, this is no longer true when the bias profile is discontinuous. In the case of a 
(5-shaped bias, Eq.([6]) leads to a discontinuity of the phase gradient in the originjj], [5]: 
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0a;(O+) — 0a;(O— ) = 7^. Exploiting again the symmetry property of the phase gradient, we 
have: 0x(O— ) =75/2 and 

0x(O+) = -7<5/2. (8) 

To summarize, an overlap junction with ABC obeys Eq.([2]) with the condition at x = //2 as 
in Eq.Q: 

0,(//2) = ha, (9) 

and the condition in the origin given by Eq.([7| or ([S]) in the cases of continuous or (5-shaped 
bias, respectively. 



SMALL JTJS 



Since the phase profile of a junction with ABC has to be even, then, for small junctions 
(/ < 1), we can use the trial function ^(x) = ax'^ + b\x\ + 00, with the parameters a and b 
to be determined from the boundary conditions and 0o treated as an integration constant. 
We will consider separately the cases of uniform bias and that of a (5-shaped bias profile [9j. 
In the former case, in order to fulfill the condition in (j7|, then b — and, to satisfy ([9]), 
a — ha/ 1, so that: (f){x) — hax'^/l + 0o- The Josephson current ij can be computed from the 
above quadratic expression: 




cos 00 



CI A / ^ 1 sin (/^o 



(10) 



in which we have introduced Fresnel's integrals defined by: sin ax^dx — 

yjTi/2aS {yj2a/Tix\ and cos ax'^dx — y/T^j2aC [ yj2a/'K x\ {a > 0). The critical cur- 



rent ic can be found by maximizing (10) with respect to 0o- Introducing the quantity 



he — yJ\ha\l/2Ti^ we get a Fresnel (or near-field) magnetic diffraction pattern: 



ic{he) 



^S\he)^C\K) 



with (j)o{he) — tan ^ C{he) / S{he). In the limit h^ 
that Zc(0) = 1. In the opposite limit, i.e., for he. 



► 0, S{he) 
00, S{he) 



(11) 

and C{he) ~ he so 
C{he) ^ 1/2, so that 
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iciha) ~ yjn jhal. Fig. 2 shows the dependence of the critical current i^, on the product hal 
(that is independent on Aj). 

In the case of (5-shaped bias, then, in order to satisfy the conditions in ([S]) and ([9]), 6 = —75/2 
and 2a = (2/ia + 7^)/^ — ij, so that: (j){x) — ijx'^/2 — 7^1x1/2 + (j)Q. When the modulus of 
the ratio lij/^s is much smaller than unity, the quadratic term can be disregarded, ij — 
means that 75 = —'^ha, i.e., the solution is piecewise linear: 



12] 



Eq.(12) works well either for small junctions (/ < 1) or near the pattern minima {ij ~ 0) 



or for large field values {h » 2). The Josephson current corresponding to Eq.(12) can be 
computed: 



cos 00 I 1 — COS 



sm 00 sm 



(13) 



The critical current ic{ha) can be found by maximizing (13) with respect to 0o, to give: 



ic{ha) = 



sin /ia//4 



(14) 



The integration constant is the sawtooth function: 



't'oiha) — tan~^ cotan 



(15) 



In other words, for a small (5-biased JTJ with ABC we expect a Fraunhofer-like magnetic 
diffraction pattern with a field-periodicity twice larger than that of a small uniformly biased 
JTJ in a uniform magnetic field /i", for which ic{h'^) — sin(/i^//2)/(/i'"//2). 



LONG JTJS 

Eq.(|2|) with ABC has been numerically integrated in order to find the magnetic diffraction 
patterns ic{ha) of JTJs having lengths larger that Xj. The details of the numerical technique 
can be found in Ref. [4^ together with the analogous results in the case of symmetric boundary 
conditions. Figs.3(a)-(d) show the numerically obtained ic{ha) for JTJs with ABC having 
normalized lengths / = 2, 4, 8 and 16, respectively: the full dots refer to the uniform bias, 
while the open dots correspond to the point injected current. Note the very good agreement 
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FIG. 2: Fresnel-like magnetic diffraction pattern ic{\ha\l) as in Eq.(ll)in which — ^/\ha\l/27T 



between the theoretical dependence of Fig. 2 and the result of the numerical simulations for 
/ — 2 shown by the dashed Hne of Fig. 3 (a). 

The validity of Eq.([5]) has been numerically checked both in the case of uniform bias and 
(5-shaped bias. In the latter case, for / = 2, the phase profile could be well approximated by 



= haX^ for |/ia| < 1 (0 < 00 < 7r/2) and by Eq.(I2) for > I [with 0o given by 
Eq.(I5)]. For / > 4 the linear approximation is good for \ha\ > 2 — 3. Generally speaking, 
it is observed that the magnetic diffraction patterns become more and more asymmetric as 
the junction normalized length increases. Parenthetically, we note that what is measured in 
the experiments is the maximum bias current that by virtue of Eq.Q can be quite different 
from the critical current. This might explain why the effects of a non-perfectly uniform 
magnetic field have never been reported in the literature so far. 

An asymmetric magnetic field drastically modifies also the dynamics of a planar JTJ. It is 
easy to see that the Josephson phase parity does not allow the standing wave resonances 
leading to current singularities such as Fiske and flux flow steps observed in presence of a 
uniform external fieldflO]. In fact, the magnetic fleld with opposite directions at the junction 
ends forces a chain of fluxons entering on one side and a chain of antifluxons entering on 
the other side, with the total magnetic flux in the barrier being always null. Furthermore, 
a small amplitude asymmetric fleld provides an asymmetric tuning of the average speed of 
one (or more) soliton shuttling back and forth along the junction, opposite to the symmetric 
tuning typical of a uniform fleld [TT]. 
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(c) (d) 

FIG. 3: Numerically computed magnetic diffraction pattern for intermediate length and long 
Josephson tunnel junctions with asymmetric boundary conditions: ic vs. ha for / = 2,4, 8, 16. 
The full dots refer to uniformly biased junctions, while the open circles correspond to 5-biased 
ones. 

CONCLUDING REMARKS 



In summary, it has been discussed how the magnetic properties of a planar Josephson tunnel 
junction change when the textbook assumption of a perfectly homogeneous magnetic field is 
reverted into a fully asymmetric one. In the most general case when the magnetic field has 
different amplitudes at the junction extremities, due to the system non-linear nature, the 
problem cannot be split in two subproblems with properly chosen symmetric and asymmetric 
boundary conditions, unless the conditions for linearizing the current-phase relationship 
occur [12j. The experimental verification of our findings has been planned. 
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